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SUMMARY 


An analytical method has been developed for studying the combined effects 
of rotational coupling and nonlinear aerodynamics on aircraft response for 
specified control inputs. The method involves the simultaneous solution of 
two nonlinear equations which are functions of angle of attack , roll rate, and 
control inputs. The method is applied to a number of maneuvers for a fighter- 
type aircraft. Time history responses verified the usefulness of the analysis 
for predicting a variety of response characteristics caused by interacting non- 
linear aerodynamic and inertial effects, including spin conditions. 


INTRODUCTION 

Maneuvering aircraft sometimes undergo sudden divergent motions not pre- 
dicted by linearized response analysis. The basic theory explaining this non- 
linear response was established by Phillips in reference 1 . The analysis of 
reference 1 predicted that divergentlike motions would occur for certain 
critical roll rates when the linearized analysis predicted satisfactory behav- 
ior. Numerous extensions have been made to Phillips' analysis (for example, 
refs. 2 to 5) . These studies showed that Phillips' critical roll rates were 
related to pseudosteady solutions of the equations of motion and that the 
divergences ("jump" phenomena) observed in simulated maneuver calculations were 
transitions from one pseudosteady solution to another. These studies consid- 
ered roll rate as the basic independent variable, and their main objective was 
to predict peak values of the incidence angles, 3 and a, in the vicinity of 
the critical values of roll rate. However, the pseudosteady values were found 
to be poor estimates of these critical peak values. The reason for this dis- 
agreement was shown in reference 6. An elegant perturbation expansion of the 
weight-component variations was used to show that these variations were very 
important in estimating peak values of a and 3 for constant roll rates near 
the critical values. The fact that weight- component variations are ignored in 
pseudosteady-state analysis explained the failure of the previous studies. How- 
ever, the method of reference 6 is also of limited usefulness in predicting peak 
values in critical rotational maneuvers because it is based on the assumption of 
constant roll rate. Many simulator studies have shown that for critical combi- 
nations of control inputs, all the coupled responses, including roll rate, are 
very irregular, so that calculations based on constant roll rate are not 
realistic. 

In reference 7, Schy and Hannah took a different approach to the applica- 
tion of pseudosteady-state (PSS) analysis. They found that the quantitative 
calculations of peak values can best be performed by computer integration of 
the nonlinear differential equations, but the pseudosteady-state analysis could 


be very useful in the qualitative analysis of critical rotationally coupled 
maneuvers and in obtaining approximate values for critical control combina- 
tions. Multiple pseudosteady solutions were plotted over the domain of 
allowable control inputs, and critical control combinations were identified 
as those for which there was a tendency to diverge from the normal response 
regime toward another solution regime corresponding to some highly coupled, 
nonlinear response. Stability characteristics of each solution were also 
shown to be useful in interpreting the pseudosteady results. Although it was 
shown in reference 7 that the method could be extended to include nonlinear 
aerodynamic effects, the results presented were based on the assumption of 
linear aerodynamics . 

The analysis of reference 7 was extended in reference 8 to include aero- 
dynamics which were nonlinear in angle of attack. Although it was assumed in 
reference 8 that the PSS solutions were mainly useful in considering rapid 
rotational maneuvers as nonlinear perturbations from horizontal flight, it was 
shown that the PSS analysis also yielded high-angle-of-attack, spinlike solu- 
tions. However, the equations of motion and aerodynamics used in reference 8 
were inadequate to represent the high-angle-of-attack regime. 

In references 9 and 10, Mehra and his co-workers have shown that catas- 
trophe and bifurcation theory can be applied to pseudosteady-state analysis of 
airplane maneuvers. They have developed a sophisticated computer program, 
based on continuation theory, for calculating both pseudosteady and truly steady 
solutions and have shown results for several example cases. The method of ref- 
erences 9 and 10 allows for nonlinear aerodynamics and is highly mathematical in 
nature. 

In the present study, the emphasis is on the physical and engineering 
interpretation of the theory and results, with a minimum of mathematical sophis- 
tication. The objective is to extend the preliminary analysis of reference 8 by 
giving a more detailed presentation and description of the results. The method 
of reference 8 is rederived using a more realistic dynamic and aerodynamic 
representation. The present study confirms that the spinlike PSS solutions at 
high angles of attack do approximate the true equilibrium spin conditions for 
the aircraft when the aerodynamics are represented realistically. Therefore, 
to the extent that variation with angle of attack adequately represents the non- 
linear aerodynamics, the PSS method provides an effective means for surveying 
maneuver characteristics and their stability properties (for any combination of 
control inputs) over the whole regime of angle of attack. If more variables are 
needed to define the nonlinear aerodynamics, then the PSS solution would have to 
be extended to include these as independent variables. 

PSS solutions are given for a variety of pitch and roll maneuvers of a 
fighter-type aircraft, and time history responses are presented to validate the 
PSS prediction procedure. 
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SYMBOLS 


All aerodynamic data and flight motions are referenced to the principal 
axis system. 


[A] 

A 0' 

a 3 * 


characteristic stability matrix, equation (12) 


Aq#A-],A 21 coefficients of polynomial equation (11) 
\ / A4 , A5 J 


B 0' B 1 \ 
B 2' B 3 J 

b 


C l 
C- 


m 


n 3 , dyn 

C X 

C Y 

c z 

c 

D 

d 0' d 2 

E a 

g 


coefficients in equation (10) defining 3 


wing span 

rolling- moment coefficient 
pitching-moment coefficient 
yawing-moment coefficient 

“ c ng ^ c 2g a 

X 

longitudinal-force coefficient 
side-force coefficient 
vertical-force coefficient 
mean aerodynamic chord 
drag 

coefficients of equation (7) 
aileron effectiveness parameter, 
acceleration due to gravity 
I x ,I Y ,I z ,I X z body-axis moments and product of inertia about center of mass 
j x = d z - i Y )/i x 

j y = d z - i X )/i y 


' n B Cz 6 Cn 6 a / Cz 6. 
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Jz 

= (Iy - I X )/I Z 


m 

mass of aircraft 


P/q f r 

body-axis rolling, pitching. 

and yawing angular rates 

Qo»Ql rQ2 

coefficients in equation (8) 

defining q 

%o 

dynamic pressure, -pv z 

2 


R 0» R 1 

coefficients in equation (9) 

defining r 


S wing area 

t time 

V airspeed 

W weight of aircraft 

a angle of attack 

0 angle of sideslip 

A determinant of matrix, equation (7) 

Y flight-path angle 

$a aileron deflection, positive when right aileron trailing edge is down 

0 e elevator deflection, positive when trailing edge is down 

A6 e elevator increment from one-g trim condition (A6 e = 6 e + 3.1°) 

<$ r rudder deflection, positive when trailing edge is left 

p air density 

4> , 0 Euler angles 

to angular velocity, (p2 + q2 + r 2jl/2 

Coefficients and derivatives: 


dci 

CZ<S a ” isl 


<3ooSb 




c l, 


l6 a 

iSa ' 


4 


I 


9C Z 

C *«r = zT t 


l & r 


%o Sh 


C ^r 




' 6r = 


9C Z 

~ aF 


Q Sb 

' 17 c ’s 


- Zg 

Zg = - 

J x 


c z„ = 


dCi 


P pb 
8 2V 


%> sh b 

Zp Cl 

P I x 2V f P 


~ z D 

1 >-Z 


9C Z 

C l = 

r rb 


9 2V 


^o Sb b 

2 r = C Z 

Iy 2V ‘c 


A Z 

Z r = ~ 
J X 


z 0 - z 6 a 5 a + 2 6 r 6 r 


< 3 0 o Sc 

m = C 


m 


m 


m = 


ac 


Cm 6 Q " 


m 


e 36. 


‘W 30 


m<S e = i y Cm <S e 


m 6 . 


m 6 ~ = 


9C 


Crn q 


m 


qc 

9 ^ 


<3co Sc c 

^ = F7 ^ Cm q 


m q- = 
q J Y 


A A A 

mg = m + mg <S € 


' n 6 


9C 
a 96- 


n 


n 6 = = 


%o Sh 

17 C " S a 


n 6- 

/s C 

n<Sa " F~ 
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c y 6 


r 


9Cy 

d&~ r 


q^S 

Y( ^r _ mV Cy 6 r 


c Yg = 


3Cy 

aT 



C Y$ 




%o S b 

Cy 

mV 2V P 


c Y r 



Y r 


b 

Cy 

mV 2V r 


Y 0 = Y 6 a 5 a + Y 6 r 6 r 


Subscript: 

a partial derivative with respect to angle of attack 

A dot over a symbol indicates a derivative with respect to time. A bar over a 
symbol indicates a pseudosteady-state solution. 

PROBLEM FORMULATION AND SOLUTION 


An analytical method for finding approximate steady-state solutions to the 
nonlinear equations of motion of an aircraft in maneuvering flight is discussed 
in this section. The assumptions used in the analysis and the equations for the 
computation of the stability characteristics of the approximate solutions are 
also discussed. 


Assumptions 

The nonlinear equations of motion used in the calculation of time history 
responses are the same as those used in reference 11 except that they are 
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written with respect to principal axes and it is assumed that speed is constant 
and that sideslip angle is small. Air density is assumed to be constant. The 
equations are given in appendix A. 

To obtain equilibrium, or steady-state, solutions, the equations of appen- 
dix A must be solved with all the time derivatives set to zero. This approach 
was used in reference 11 to calculate equilibrium spin conditions for air- 
craft. (The computer program used in the study of ref. 11 is described in 
ref. 12.) The resulting motions calculated in reference 11 consisted of heli- 
cal paths about a vertical axis. Pseudosteady states (PSS) differ from these 
truly steady states in that Euler-angle variations which determine the weight- 
component variations with respect to the body axes are ignored. This approxima- 
tion is appropriate in analyzing rotational coupling effects in rapid maneuvers 
because it removes the constraint that the motion be about the vertical axis. 
This was achieved in the present analysis by ignoring the 6 and <j> equations 
of appendix A and by assuming that 0 = a (Y = 0) and 0=0 in the a and 
3 equations. The PSS analysis also ignores the V equation of reference 11 by 
assuming that a constant equilibrium cruise velocity is maintained. (Modifica- 
tions to the assumptions of constant velocity, 0 = a, and 0=0 are made in 
a later section of the report which gives comparisons between PSS solutions and 
true equilibrium spin conditions.) The resulting five equations, with their 
time derivatives set to zero, are used in the following analysis to obtain PSS 
solutions for any combination of control inputs. The basic assumptions of the 
PSS method of analysis are (1) that important qualitative dynamic characteris- 
tics of rapid aircraft maneuvers can be obtained by studying the state-space 
loci and stability properties of the approximate equilibrium solutions over the 
range of practical control combinations and ( 2 ) that critical control combina- 
tions can be identified by significant changes in the patterns of the state- 
space loci. 


Solution for Pseudosteady States 

Equations (A1 ) to (A5) of appendix A could be solved for steady states 
using an optimization algorithm such as that used in reference 11 to compute 
equilibrium spin conditions. The method of reference 11 required a five- 
parameter search. The simpler PSS formulation permits reduction of the problem 
to a two-parameter search, so that it is possible to develop the following 
simple and efficient PSS solution procedure. 

Equations (Al) to (A5) can be written as 


pq - Bng - rn r = n 0 + pn p 

( 1 ) 

qniq + pr = -m 0 

( 2 ) 

r(cos a - Y r ) - 8 Y 3 = Y 0 + (Y p + sin a)p 

(3) 
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3 2 g + P^p r ^ r ^ 0 


(4) 


q - 3 (p cos a + r sin a) + Zq + g/V = 0 (5) 

The overbars in the preceding equations indicate PSS solutions and, for conve- 
nience, the arguments for the aerodynamic coefficients have been dropped. 

Equations (1), (2), and (3) are linear in q, 3, and r and can be solved 
for these variables as functions of p, a, and control inputs. These equations 
can be written in matrix form as follows: 


~~ 

/V 

A 






A 

p 

” n r 

" n B 

q 


n 0 


"P 

A 

m q 

P 

0 

r 

= 

1 

B> 

0 

+ p 

0 

0 

cos a - Y r 

- y b 

3 


K 

0 

l 


Yp + sin a 


The determinant of the preceding matrix is given by 


A = D 2 p 2 + D 0 


(7) 


where 


D 2 = -Yp 


Dq = -fhg[ Ygn r + ng(cos a - Y r )] 


Solutions for q, r, and 3 are as follows: 


q = (Q2P 2 + QiP + Qq)/^ 


r - (R-|P + 'Rq) 


B = (B^p^ + B 2 P 2 + B-jp + Bq)/A 


( 8 ) 

(9) 

( 10 ) 
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where 


Q 2 = ng (Yp + sin a) - Ygn p 
Ql = ngY 0 - n 0 Yg 


Q 0 = m 0 [Ygn r + ng(cos a - Y r ) ] 

R] = m 0 Yg - m q [ng (Y p + sin a) - Ygn p ] 

R 0 = i"q(n 0 Yg - ngY 0 ) 

B3 = Yp + sin a 
B 2 = Y 0 

B-] = (cos a - Y r ) (mg - mqn p ) + Siqn r (Y p + sin a) 

Bg = mq[n r Yo + ng(cos a - Y r ) ] 

Substituting these values for q, r, and 6 into equation (4) gives the fol- 
lowing polynomial in p: 

A 3 p5 + A 4 p^ + A 3 p3 + A 2 p 2 + A-jp + Ag = 0 (11) 

where 

A 5 = D 2 (2gB 3 + i p D 2 ) 

A 4 = D 2 (ZgB 2 + ZgD 2 ) 

A 3 = *g( D 0 B 3 + D 2 B] ) + D 2 (Z r Ri + 2^ p Dg) - Q 2 R-| 

a 2 = ^3 ( d 0 b 2 + d 2 b 0^ + R 0^r D 2 “ Q2) + 2 ^0 D 0 D 2 - Ql R 1 
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A-] = Do(ZgB-) - l r R-| ) “ QiRo ~ Qo R l + ^pDo 2 

A A A 

A 0 “ Dq(^3 Bq + Z r Ro +z O D o) " Qo R o 


Simultaneous solution of equations (5) and (11) yields PSS values for p, 
a, q, r, and 3 for specified values of the control inputs. (Eqs. (8) , 

(9) , and (10) are used to compute q, r, and 3.) Various iterative search 
algorithms could be used to solve equations (5) and (11). The search procedure 
employed in the present study is described in appendix B. 

Stability of PSS Solutions 

The stability of any particular set of PSS solutions was calculated by 
linearizing the first five equations of appendix A about the PSS solution. The 
matrix whose characteristic roots determine the stability of a perturbed motion 
is given by 


where 


n r 

- J ZP 


n B 

n p " 


JyP 

H'q 


0 

JylT 

<3a 

- cos a 

0 


*3 

Yp + sin a 

3a 

- J xi 

-J x r 


*3 

Z P 

Pa 

sin a 

1 

-p 

cos a -r sin 

a -8 cos a 

a a _ 


= n 3a e + 

w + 

n r a " + n °a 



• 

<3a 

= ^ + 

m °a 




3 a 

= r sin a + p cos a + Yg B + 

Y Pa^ + Y r a r + Y °a 


• 

Pa 

- ! 6a B + 

2 Pa? + 

Z r r + l n 
r a °a 



“a 

= p3 sin 

a - rB 

cos a + Zn 
u a 




= [a] (12) 
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In the preceding equations, the absence of carets (~) on the moment derivatives 
indicates that the inertia ratios are not factored out as in the equations of 
appendix A. 

It should be noted that the stability characteristics as defined by [a] 
are approximate, since the varying weight components (in the a and 3 equa- 
tions of appendix A) are neglected. Also, stability derivatives in column 5 
of [a] were computed by taking numerical slopes using linear interpolation on 
the nonlinear aerodynamic data. This procedure introduces discontinuities in 
the stability derivatives at breakpoints in the tabular aerodynamic data used 
in the computer program. However, these approximate stability characteristics 
are considered adequate for the present analysis since the PSS solutions are 
only approximate steady states. 


DESCRIPTION OF AIRCRAFT CONFIGURATION 


The assumed aircraft of the present study is representative of a twin- 
jet, swept-wing fighter. The aircraft was assumed to be flying at an altitude 
of 13 720 m and at a Mach number of 0.9 (V = 266 m/sec, q^ = 8400 Pa). It 
should be noted that this flight condition is more representative of a flight 
test condition than of normal operational flight. 


Weight, inertia, and dimensional characteristics for the aircraft are 
given in table I. Aerodynamic data used in the analysis were taken from refer- 
ences 13 and 14 and are tabulated in table II. With the exception of C 
C n g, and c yg/ the 3=0 data of reference 13 were used to represent the 
coefficients. The coefficients C c ng' and c Yg were determined at each 
angle of attack by taking the slope P of the data of reference 13 between side- 
slip angles of ±5°. (It should be noted that the assumption of linearity in 
3 can introduce errors in time history responses which involve large transient 
motions in sideslip angle.) The dynamic stability derivatives were taken from 
reference 14. Linear interpolation was used between the data points given in 
table II. Variations of several important lateral-directional parameters with 
a are given on figure 1. These include c ng' C ?g' and two P a rameters which 


are strongly dependent on these 


Cn 3 dyn' wh ^ clri -*- s commonl y used as an indica- 


tor of divergence tendency, and E a , which is an approximate measure of aileron 
roll effectiveness. 


RESULTS AND DISCUSSION 

The pseudosteady-state (PSS) method has been applied to a number of maneu- 
vers for the assumed aircraft. These included aileron roll maneuvers about a 
one-g condition (Lift = Weight, Pitching moment = 0) and roll maneuvers initi- 
ated from pitch-up (Acceleration > g) and pitch-down (Acceleration < g) condi- 
tions. The trim conditions considered in the analysis are summarized in 
table III. PSS predictions for the aforementioned maneuvers are now discussed. 
Time history responses are then presented to validate the PSS predictions. 
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PSS Solutions 


PSS solutions for angle of attack, roll rate, yaw rate, sideslip angle, and 
pitch rate for various aileron inputs are shown on figure 2 for the one-g trim 
condition. (The one-g condition actually required S e = -3.1°, which is here- 
after referred to as AS e = 0°; all other elevator settings are given as incre- 
ments from the one-g condition.) Figure 2(a) shows angle-of-attack solutions. 
The curves are numbered from 1 to 6 in order of increasing angle of attack. The 
curves are also labeled as P or N to indicate positive or negative roll rates, 
as shown on figure 2(b). The stability of the solutions is indicated on the 
figure; the stability characteristics are tabulated in table IV for selected 
aileron inputs. Consider the <5 a = -25° results of table IV for the 2P and 3P 
curves. Figure 2 shows that the values for p, a, (3/ q, and r ace approx- 
imately the same at this condition. Therefore, we would expect the character- 
istic roots to be approximately the same for 6 a = -25° on the 2P and 3P 
curves. The differences shown in table IV result from the previously mentioned 
method for finding the stability derivatives in column 5 of [a] since 2P and 
3P are on opposite sides of a breakpoint in the tabular aerodynamic data at 
a = 35°. (Column 5 of [a] was computed by taking numerical slopes using linear 
interpolation between data points. This procedure introduces discontinuities 
at breakpoints in the tabular aerodynamic data.) If a more accurate stability 
representation is desired, a higher order interpolation scheme should be used 
in fitting the aerodynamic data. 

Figure 2 shows that there are multiple PSS solutions for all aileron set- 
tings. As shown in previous studies (refs. 7 and 8) , the existence of multiple 
PSS solutions can lead in various ways to predictions of a divergence in 
response at certain control input combinations. Divergences are associated with 
the behavior of the basic solution, that being the normal aircraft response 
curve (curve 1, fig. 2) starting at p = 0 and <5 a = 0. If the basic solution 
disappears at some control combinations by joining another solution (curves 
2 and 3 join on fig. 2) , then the response would be expected to diverge to the 
vicinity of another PSS solution. Also, if the basic solution comes close in 
state space to another solution, a perturbed response may enter the domain of 
attraction of the other solution and a divergence might occur. 

The stability of the PSS solutions is important in determining possible 
attractor states. Solid lines on figure 2 represent stable solutions which are 
strong attractors. Short dashed lines represent solutions for which a diver- 
gence would occur. Short-long dashed lines are potential attractors since they 
indicate an unstable oscillation. Of course, the regime of each PSS solution 
has a significant effect on the aircraft response when a maneuver enters that 
regime depending on the stability properties in its domain of influence. For 
example, if the maneuver enters the domain of influence of a strongly divergent 
solution, the response should rapidly diverge to a more stable regime. 

We now analyze figure 2 with respect to possible divergences in aircraft 
response. Since the basic solution for AS e = 0 remains stable for all aile- 
ron deflections, the only possibility of a divergence would be for the basic 
solution curve to enter the domain of attraction of other curves. This occurs 
for large aileron inputs in the roll-rate solution for curves 4P, 5P, and 6P 
(fig. 2(b)). However, the basic solution is far away from the other 4P, 5P, 
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and 6P solutions of figure 2, particularly the angle-of-attack solutions. Since 
angle of attack is such an important parameter in determining the response char- 
acteristics with nonlinear aerodynamics, a divergence would be unlikely; that 
is, the aircraft response would be expected to attain the levels of the basic 
solution for all aileron inputs. Nevertheless, there is enough uncertainty in 
this case that the conclusion is verified with time history responses. 

The high-angle-of-attack solutions (a > 25°) of figure 2 are of interest 
since they resemble possible spin conditions (a above stall, large rotation 
rates). These solutions, which are typical of those obtained with other ele- 
vator inputs, are discussed in some detail in a later section of the report. 
However, for the present, only the PSS solutions in the normal operational enve- 
lope (a < 25°) of the aircraft are considered. 

Figure 3 shows PSS solutions for roll rate, sideslip angle, angle of 
attack, pitch rate, and yaw rate as a function of aileron deflection for various 
elevator increments above and below the one-g trim value (A6 e = 0°) . Only the 
basic solution curves and those associated with the basic solution are given on 
figure 3. No solutions for angles of attack greater than 25° are shown. The 
stability characteristics at various points on the curves of figure 3 are given 
in table V. 

Consider the one-g (AS e = 0) solution. Figure 3(a) shows. that as aileron 
deflection is increased, roll rate begins to saturate; that is, the increasing 
positive sideslip angle (fig. 3(b)) acts through C to reduce roll rate. 

This effect is even more evident for the pitch-up maneuvers with A6 e = -4° 
and -8°. At the high angle of attack corresponding to A6 e = -8° (fig. 3(c)), 
the aileron effectiveness in producing roll almost vanishes and at A6 e = -12°, 
a roll reversal occurs (negative aileron deflection generally gives positive 
roll rate) . This is predicted by the aileron effectiveness results of fig- 
ure 1 which show a reversal in effectiveness for angles of attack greater than 
about 20°. 

Now consider the pitch-down (A6 e > 0°) rolling maneuvers shown on figure 3. 
For A6 e » 6° and 8°, there are multiple PSS values for certain aileron inputs. 
For these cases, as aileron deflection increases beyond a certain level, the 
basic solution disappears by joining a divergent PSS solution. For aileron 
inputs greater than the value at which the basic solution disappears, the air- 
craft response would be expected to diverge to the vicinity of other stable 
PSS solutions. 

Consider the PSS solutions for A6 e = 8° shown on figure 3. This case 
corresponds to rolling maneuvers about a -0.6g condition. For aileron inputs 
up to about -9°, the steady-state response would be expected to reach the 
levels of the basic solution. As the aileron deflection is increased from -9° 
to about -10°, the following would be expected to happen. The roll rate 
(fig. 3(a)) would jump from about 70 deg/sec to about 145 deg/sec; angle of 
attack (fig. 3(c)), from about -6° to about 7.5°; and yaw rate (fig. 3(e)), 
from about -8° deg/sec to about 20 deg/sec. The steady-state values for side- 
slip angle and pitch rate (figs. 3(b) and 3(d)) would be about the same as 
before the divergence. Similar predictions could be made for AS e = 6° on 
figure 3. 
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The roll-rate response predictions of figure 3(a) for A6 e = 4° are inter- 
esting. This case, which corresponds to rolls about a nearly zero-g condition, 
is the approximate boundary between different types of coupled responses. Fig- 
ure 3(a) shows that for aileron deflections up to about -26°, the roll rate 
appears to be saturating. However, for even small increases in aileron deflec- 
tion above -26°, the roll response increases rapidly. This roll response was 
predicted by the analysis of reference 1 . That analysis predicted an unstable 
aircraft response if the roll-rate magnitude was between the values for the 
undamped natural frequencies of the longitudinal and lateral modes of the air- 
craft.^ The longitudinal and lateral frequencies are approximately (-mg^) 1 ^ 
and (ng)V2 f respectively. For A6 e = 4°, and 6 a = -26° on figure 3, the 
longitudinal and lateral frequencies are approximately equal, both being about 
125 deg/sec. 

In contrast to the pitch-down elevator responses just discussed, consider 
the other extreme, AS e = -12° on figure 3. This case corresponds to rolls 
about a 2.2g pitch-up condition. For small aileron inputs, there is a reversal 
in aileron effectiveness and the basic solution shows an unstable oscillation. 
For aileron deflections greater than about -1.5°, the basic solution disappears 
and figure 3 shows that there are no other solutions for the angle-of-attack 
range shown. For all aileron inputs, there are, however, several high-angle- 
of-attack solutions similar to those shown on figure 2. Thus, even for small 
aileron inputs, we would expect the A 6 e = -12° response to diverge toward the 
high-angle-of-attack PSS solutions. 


Summary of PSS Solutions 

Three types of maneuvers were analyzed in the previous discussion. These 
included rolls about a one-g trim condition and rolling maneuvers using pitch- 
down and pitch-up elevator deflections. Figure 4 illustrates the general nature 
of the PSS solutions for these three types of maneuvers and the possible prob- 
lems associated with each. The circled numbers represent the different solution 
types, ordered by increasing angle of attack, as on figure 2. Shown on fig- 
ure 4(a) are all positive roll-rate solutions for the combinations of A6 e and 
6 a given on figure 3. Showing the solutions as a function of angle of attack 
gives a better understanding of the relationship between various solutions. The 
stability characteristics of the solutions are also shown. Negative roll-rate 
solutions would be a mirror image of the positive solutions shown on figure 4(a); 
that is, the positive curves of figure 4(a) could be rotated about the p = 0 
axis to become negative solutions, in which case the sign of the aileron deflec- 
tion would be reversed. 

Consider the pitch-down (AS e > 0) cases shown on figure 4(a). For 
AS e = 6° and 8°, the basic solution disappears at some critical aileron input. 

For aileron deflections greater than these critical values the response could 
be expected to diverge to the high-roll rate, low-angle-of-attack solutions in 
the upper left-hand part of figure 4. This region is a spin stabilized region, 
where gyroscopic momentum adds to the "stiffness" in a and 3. There is a 
possibility that the response would be attracted to the high-angle-of-attack 
solutions once the basic solution vanishes. This possibility is investigated 
with time history calculations in a later section. 
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Now consider the basic solution (1) curves of figure 4(a) for pitch-up 
maneuvers (AS e < 0) . The previously shown roll-rate saturation with increasing 
aileron deflection is evident. Also evident is the reduced aileron effectiveness 
as pitch-up elevator deflection (-A6 e ) is increased. That is, as the trim angle 
of attack corresponding to pitch-up elevator deflection increases, the roll 
effectiveness of the aileron decreases. As shown on figure 4(a), aileron effec- 
tiveness becomes zero and then reverses for angles of attack above about 20°. 

The reversal is not caused by a change in sign of , but results from a 

3 


change in sign of the aileron effectiveness parameter shown on figure 1 . The 
roll reversal is clearly seen for A6 e = -12° on figure 4(a). For this case 
the basic solution joins a divergent solution for aileron deflections greater 
than about 1.5° in magnitude and, for larger aileron deflections, the only 
attractors are solutions 3, 4, 5, and 6, which have angles of attack greater 
than about 25°. This suggests that for this case, the aircraft response would 
diverge to high angles of attack for aileron inputs greater than about 1.5° in 
magnitude. 


The high-angle-of-attack, spinlike regimes are better illustrated by show- 
ing total angular velocity as a function of angle of attack. This is shown on 

figure 4(b) for the range of A6 e and 6 a given on figure 4(a). For the 

curves shown on figure 4(b), angular velocity is composed almost entirely of 
p and r with the q contribution being insignificant. 

A comparison of figure 4(b) with 4(a) shows that for angles of attack up 
to about 25°, angular velocity is dominated by roll rate. For a > 25°, there 
is a gradual increase in yaw rate, and at about a = 45°, the p and r con- 
tributions to angular velocity are about equal. For higher angles of attack, 

yaw rate is the dominant term in the angular velocity. The 5- and 6- type solu- 
tions shown on figure 4(b) are composed almost entirely of yaw rate. 


The divergent 2-type curves shown on figure 4 may represent the spin depar- 
ture region since the aircraft response would be expected to diverge rapidly 
through this region. The oscillatory 3- type and divergent 4- type curves resem- 
ble steep and intermediate spin conditions. The oscillatory 5- type and diver- 
gent 6-type solutions may represent flat spin conditions. We refer to these as 
spin conditions, although if the aerodynamic representation were accurate, one 
would not expect the aircraft to spin at a condition which is predicted to be 
divergent. It should also be noted that assuming the velocity to be the con- 
stant cruise value is invalid at high angles of attack. This is discussed in 
a later section of the report. 


In summary, figure 4 shows that depending on the control combinations, the 
aircraft response can be attracted to three distinct regions. It can remain at 
low angles of attack and moderate roll rates, it can diverge to low angles of 
attack and high roll rates, or it may diverge to high angles of attack and large 
rotation rates. 
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Time History Responses 

Time history responses are now presented to validate the PSS predictions 
of figures 3 and 4. Responses are given for maneuvers about a one-g trim condi- 
tion, pitch-down conditions, and pitch-up conditions. While the PSS analyis 
neglected the varying gravity components, the time histories were calculated 
using the complete equations of motion given in appendix A, unless otherwise 
noted. 

One-g t rim maneuvers .- Time history calculations for rolling maneuvers 
about the one-g trim condition ( AS e = 0) agreed well with the PSS predictions of 
figure 3. At any aileron deflection, the steady-state value for all variables 
was close to that predicted by the basic solution of figure 3. Roll-rate and 
sideslip-angle responses are shown on figure 5 for two aileron inputs. PSS pre- 
dictions are also noted. As can be seen, the responses are close to the PSS 
predictions. Similar agreement was found for the other variables. One of the 
roots predicted by the stability characteristics of table V is evident in the 
time histories of figure 5. A mode having a root with a period of 2.8 seconds 
was predicted for = -5 Q . The corresponding p and 3 histories have a 
period of about 3 seconds. A mode having a root with a period of 2.1 seconds 
was predicted for <$ a = -30°. The corresponding p and 6 histories contain 
an oscillation with a period of about 2 seconds. 

Pitch- down maneuvers .- Time history responses for pitch-down elevator 
deflections (A<5 0 > 0°) were in agreement with the PSS predictions of figures 3 
and 4. For example, for A6 0 = 8°, figure 3 predicts a divergence for an aile- 
ron input of about -9.5°. Time history responses showed that the divergence 
actually occurs for aileron deflections of -10° or greater. 

Figure 6 illustrates the aircraft response before and after a coupled 
divergence for the A6 0 = 8° maneuver. For this case, which is typical of cou- 
pled responses, aileron deflection was maintained at -8° for the first 5 seconds 
and then increased to -12°. PSS predictions for each aileron deflection are 
also indicated on figure 6. The roll-rate, pitch-rate, yaw-rate, angle-of- 
attack, and sideslip-angle responses at 5 seconds are in general agreement with 
the PSS basic solution of figure 3 for S a = -8°. When the aileron deflection 
is increased to -12° after 5 seconds, the tendency to diverge is evident first 
in p and 3. The increasing negative sideslip acts through the dihedral 
effect C2g to accelerate the rolling. At about 13 seconds the coupling terms 
in the r and ot equations (appendix A) , that is the -pq and -p3 cos a 
terms, cause r and a to diverge toward the positive values predicted on fig- 
ures 3(e) and 3(c). At about 16 seconds, a similar coupled divergence occurs in 
pitch rate (caused by .the pr term in the q equation) . An interesting fea- 
ture of the responses shown on figure 6 is that all the variables initially 
diverge toward values consistent with the divergent autorotational values 
(A6 0 = 8°, <$ a = 0°) shown on figure 3. It appears that the transition between 

the two stable conditions is characterized by a passage through the divergent 
A<$ e = 8° and <$ a = 0° domain. 
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The major factors contributing to the effects shown on figure 6 are given 
on figure 7. Shown are the dominant terms in the p, r, and 3 equations. 

The p components given on figure 7(a) show that the divergence to negative 
sideslip angles shown on figure 6 is caused by the p sin a coupling term. 
Consider the components of the roll acceleration shown on figure 7(b). As can 
be seen, Zg in combination with the applied rolling moment (Zg), overcomes 
the damping in roll (Z p ) to increase the roll acceleration. The qr coupling 

term, which is not shown has an insignificant effect on p. In contrast, the 

coupling term -pq in the yaw-acceleration response of figure 7(c) is the driving 
term. Similar roll coupled effects occur in the q and & equations because 
of pr and p3 cos a. 

While the initial rapid divergence associated with rotational coupling is 
the important aspect of the present analysis, it is of interest to check the PSS 
predictions following a divergence. Continuation of time histories of figure 6 
beyond 20 seconds would show the responses converging to values near the PSS 
predictions. This is illustrated better by using an aileron deflection further 
removed from the critical value of about -10°. Shown on figure 8 are p, a, 
and r responses for A 6 e = 8° and 6 a = -15°. As shown, the responses after 
divergence are close to the PSS predictions of figure 3. 

Suppose that a time history is initiated at PSS values corresponding to the 

A6 e = 8° and 5 a = -15° condition of figure 3 and that the aileron deflection 

is then reduced to zero. The results of figure 3 suggest that the aircraft 
states would not return to their trim values (table III) but would remain at 
their autorotational values. For example, roll rate might be expected not to 
return to zero but to attain the autorotational, <$ a = 0° value of about 
140 deg/sec shown on figure 3(a). The responses are attracted to the autorota- 
tional values if the aileron is ramped back to zero at a rate of about 5 deg/sec 
or less. Figure 9 shows two roll-rate responses initiated from the A6 e = 8°, 

6 a = ~15°, PSS condition for which the aileron deflection was ramped down to 
zero at a rate of 5 deg/sec. The responses were calculated both with and with- 
out the varying weight components as shown on the figure. Both responses are 
attracted to the autorotational, PSS solution. 

However, when the aileron deflection was immediately set at zero, the vari- 
ables returned to the levels of the trim condition after about 6 seconds. This 
is because of the inclusion in the time history calculations of varying weight 
components which were omitted in the PSS analysis. When the weight-component 
variations were omitted from the time history calculations (<j> = 0 = 0) and 6 a 
was set immediately to zero, the responses continued to oscillate about the 
autorotational values. 

Pitch-up maneuvers .- Time history calculations were in agreement with the 
PSS predictions of figure 3 for AS e = -4° and -8°. The steady-state values 
for all variables attained the level of the PSS predictions for all aileron 
inputs with no tendency to diverge. 

Now consider the results shown on figure 3 for AS e = -12°. The low roll 
rates associated with this maneuver are not of the type for which the analysis 
was intended. However, the maneuver is of interest since the results of fig- 
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ures 3 and 4 predicted roll reversal for the basic solution and a divergence to 
high angles of attack and large rotation rates for small aileron deflections. 

Time history responses verified the PSS predictions for A6 e = -12°. For 
all aileron deflections, the response initially diverged to the vicinity of the 
high-angle-of-attack PSS predictions. This behavior is illustrated on fig- 
ure 10. Shown are p, a, £, q, and r time histories for an aileron deflec- 
tion of -5°. This response is typical of those for other aileron deflections. 
The figure shows that during the first 3 seconds, a, 3/ and r gradually 
increase in magnitude. At a time of about 3 seconds, there is a rapid diver- 
gence in all response variables. Continuation of the time history beyond 
5 seconds shows the response approximately repeating the history shown on fig- 
ure 10 between about 3 to 5 seconds, with p and r oscillating between large 
positive and negative values. 

The initial divergence shown on figure 10 is caused by both aerodynamic and 
coupling effects and can be understood by examining the component parts of the 
equations of appendix A. The dominant components of the 3, r, and p equa- 
tions are shown on figure 11. Figure 11(a) shows that the increasing sideslip 
angle is caused by the r cos a coupling term, while figure 11(b) shows that 
the unstable ng causes yaw accelerations to diverge to negative values. That 
is, C n o becomes negative for angles of attack greater than about 22° (fig. 1) 
and acts through 3n R to increase the magnitude of yaw rate. This is reflected 
in the stability predictions of table V which show unstable oscillatory roots 
for A6 e = -12°. 

The components of the rolling-moment equation are shown on figure 11(c). 

For about the first 3 seconds r 2 r , 3Zg, and 2 q nearly balance each other 
to result in a small negative roll rate. At slightly less than 3 seconds, the 
dihedral effect Ci^ causes roll acceleration to diverge to large negative 
values. As seen on P figure 11(c), 3 Zg is positive between about 2.2 seconds 
and 2.9 seconds. Since 3 is positive during this period (fig. 10), Cjg 
must also be positive. This is, in fact, the case between angles of attack 
of about 25° and 35° as shown on figure 1 . When angle of attack exceeds 35° 
at about 3 seconds, becomes negative and causes roll acceleration to 

diverge as shown on figure 11(c). 

The responses shown on figure 10 for the first 3 seconds are largely due to 
aerodynamic effects. However, the continuing rapid oscillations after 3 seconds 
are highly coupled about all aircraft axes. 

It should be emphasized that the physically significant portion of the time 
history of figure 1 0 is the initial divergence occurring during the first 3 sec- 
onds. After that, angle of attack is above the stall and the constant velocity 
assumption is invalid since velocity would decrease rapidly at high angles of 
attack. In addition, for large sideslip angles, the linearity assumption in 
such coefficients as Cj is also invalid. However, it is of academic interest 
to examine these responses to see how they are influenced by the spinlike PSS 
solutions obtained in the analysis. 

A clearer picture of how the response of figure 10 is attracted to the 
high-angle-of-attack solutions is obtained by showing p and a, r and a. 
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and id and a phase plots for the time history of figure 10. The first 9 sec- 
onds of these trajectories is given on figure 12. Also shown are the PSS 
predictions. The predictions are numbered 2 through 6 in accordance with the 
notation of figure 4 and the stability of the predictions is also indicated. 
Consider the p and a and r and a responses shown on figures 12(a) and 
12(b). Initially the motion seems to be attracted toward PSS solutions 3 and 4. 
Then there is a rapid divergence toward solutions 5 and 6, followed by a return 
to the lower-angle-of-attack region. For the remainder of the trajectory, both 
responses are in the domain of influence of solutions 2, 3, and 4. The go and 
a response shown on figure 12(c) follows a pattern similar to these, except 
that it ini tally passses through the divergent solution 2. Continuation beyond 
9 seconds shows the responses repeating similar patterns in p and a, r and 
a, and go and a space. 


Comparison of PSS Solutions With Equilibrium Spin Conditions 

As previously shown on figures 2 and 4, PSS solutions which occur at high 
angles of attack resemble possible spin conditions in that they involve high 
rotation rates. Consider again the results of figure 4. The divergent 2- type 
curves may represent the spin departure region through which the aircraft 
response would be expected to diverge rapidly. The oscillatory 3-type solutions 
resemble steep spin conditions while the divergent 4-type solutions may repre- 
sent intermediate spin conditions. The oscillatory 5- type and divergent 6- type 
solutions resemble flat spin conditions. Comparisons are now made between PSS 
solutions and equilibrium spin solutions obtained using the method of refer- 
ence 11. Before making these comparisons, the following points should be empha- 
sized. In the method of reference 11, the motion occurs in a helical path about 
the vertical. In addition, all seven equations of appendix A, along with the 
equation for rate of change of velocity, are used in the method of reference 11. 
In brief, the method uses a minimization algorithm to find those combinations of 
trajectory variables for which all time derivatives of the equations of motion 
vanish. 

The PSS analysis neglects weight-component variation effects, so that 
pseudosteady-state solutions other than helical motions about a vertical axis 
are permitted. In a truly steady motion, such as an equilibrium spin, 0 and 
4> are fixed. Therefore it is reasonable to expect that approximate spin solu- 
tions may be found by the PSS analysis, since in the true spin the weight com- 
ponents are constant, even though changed from PSS values because 0 ^ a and 
<)> 7 * 0. Also, the moment balance in which the weight plays no role is more 
important than the force balance in a spin. 

The PSS analysis shown up to this point assumes constant velocity regard- 
less of angle of attack. In a true equilibrium spin, velocity decreases rapidly 
as angle of attack increases above the stall. A good approximation of the true 
velocity in a spin can be obtained by calculating the velocity for which the 
drag is equal to the aircraft weight. in a true equilibrium spin, 0 is not 
equal to a (y / 0) as was assumed in the PSS analysis. A reasonable 
approximation would be to assume that y = -90° (0 - a = -90°) and 4> = 3 at 

high angles of attack. 
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An additional point which should be emphasized is that the aerodynamic 
representation of this report may be unrealistic for investigating spinning con- 
ditions, Numerous studies (ref. 15, for example) have shown that a more realis- 
tic representation would be for the aerodynamic data to be functions not only of 
a and B, but also of spin rate. However, as was shown in reference 11, aero- 
dynamic data of the form used in the present study can give reasonable predic- 

tions of true equilibrium spin conditions for fighter-type aircraft. 

A comparison of high-angle-of-attack PSS solutions with equilibrium spin 
solutions is given in table VI for the one-g, 6 a = 0 condition. These 
results are typical of those obtained with other control deflections. Shown 
in table VI are the predictions using the method of reference 1 1 including the 
eighth order stability characteristics. Also shown are three sets of PSS pre- 
dictions and their fifth order stability characteristics. PSS solutions were 
first calculated by assuming that Mach number was 0.9, 0 = a, and (j) = 0. 

They were then calculated using the velocity for which the aircraft drag was 
equal to the weight (D = W) with 0 = a and § = 0. The final PSS calcula- 
tion assumed that D = W, <{) = 3/ and y = -90° (0 - a = -90°) . An examination 

of the a and B equations of appendix A shows that assuming 4> = B and 
y = -90° eliminates the g/V term in the a equation (for small B) and intro- 
duces a (g/V) (sin a - 1)3 term in the B equation. This was implemented in 
the PSS analysis by replacing Yg in equation (3) with Yg - g/V(sin a - 1) 
and setting g/V to zero in equation (5). As shown in table VI, both methods 
yield four high-angle-of-attack solutions. For the lowest and two highest 
angles of attack, there is reasonable agreement between the predicted values of 
the variables when the velocity corresponding to D = W is used. A slight 
additional improvement is obtained using D = W and y = -90°. The PSS pre- 
dictions of table VI for the intermediate spin at an angle of attack of about 
60° are not in as good agreement with significant variations occurring in all 
variables. The fifth and eighth order stability predictions shown in table VI 
follow similar trends in all cases. The results of table VI suggest that the 
high-angle-of-attack solutions obtained in the PSS analysis are reasonable 
approximations to the equilibrium spin conditions obtained using the method of 
reference 11. 


CONCLUDING REMARKS 

A pseudosteady-state (PSS) method was presented for studying the combined 
effects of rotational coupling and nonlinear aerodynamics on aircraft response. 
The method involves the simultaneous solution of two nonlinear equations which 
are functions of angle of attack, roll rate, and control inputs. Emphasis was 
placed on the physical. and engineering interpretation of the method and results, 
with a minimum of mathematical sophistication. 

The method was applied to a number of pitch and roll maneuvers for a fighter 
type aircraft. Included were aileron rolls initiated from a one-g trim condition 
pitch-down elevator conditions, and pitch- up elevator conditions. 

For pitch-down rolling maneuvers, divergence was found to occur for moder- 
ate aileron inputs. The resulting response was attracted to a stable PSS regime 
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characterized by high roll rate and low angles of attack and sideslip. For 
pitch-up rolling maneuvers the dominant effects were the degradation of lateral 
stability and control characteristics at high angles of attack. The combination 
of roll reversal and lateral divergence resulted in large, irregular, highly 
coupled responses. For example, for certain combinations of pitch-up elevator 
and aileron deflections, it was shown that there were no stable PSS solutions, 
and the resulting violent motions could be qualitatively understood as a cyclic 
transition between neighboring PSS solutions. Very high angles of attack and 
high rotation rates in this response resulted from the influence of spinlike PSS 
solutions. 

Four spinlike PSS solutions were compared with true equilibrium spin con- 
ditions. The steep and two flat spin solutions were shown to correspond well 
to the exact solutions, whereas the intermediate, divergent solution was less 
accurate. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
October 22, 1980 
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APPENDIX A 


EQUATIONS OF MOTION 

The equations of motion used in the calculation of time history responses 
are listed below. The equations are written with respect to principal axes. 

Speed and air density are assumed to be constant and the sideslip angle is assumed 
to be small. 


r = J z [-pq + 3 %(a) + p n p (a) + r n r (a) + 

q = J Y [pr + q fftq(a) + m 0 (a,6 e )] 

3 = -r cos a + p sin a + 3 Yg(a) + p Y p (a) 
g 

+ -[cos 9 sin 6 + (cos a sin 0 - sin a 
V 


n 0 ( a ' S a , 1 (AT ) 

(A2) 

+ r Y r (a) + Y 0 (a,S a ,<5 r ) 
cos 9 cos <}>) 31 (A3) 


p = J x [-qr + 3 *2 3 (ex) + p 2 p (a) + r Z r ( a) + Z 0 ( a ' 6 a' 5 r)] 

a = q - 3 (p cos a + r sin a) + Z 0 (a,6 e ) 
g 

+ -(cos a cos S cos <|) + sin a sin 0) 

<£> = p + tan 0 (q sin <)> + r cos 4>) 

0 = q cos 4> - r sin (j> 


(A4) 


(A5) 

(A6) 

(A7) 


where 


no(a.,S a ,6 J .) = n<$ (a) 6 a + n* (a) 6 r 
a r 

mo( a ' 5 e) = m(a) + n< 5 e (a) 6 e 
Y 0 (a,5 a ,6 r ) = Y<s a (a) 6 a + Ys r (a) 6 r 
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APPENDIX A 


A 

^ 0 1 ) 

= 2s a (a) <$ a + ?s r (a) 6 r 


Z 0 (a,<5 e ) = 

z (a) + z§ e (a) 6 e J cos a - 

X(a) + X(5 e (a) 6 e 

—i 
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APPENDIX B 


SOLUTION PROCEDURE 

In order to find PSS solutions for p, a, q, r, and 3, it is necessary 
to obtain a simultaneous solution to equations (5) and (11), given in the sec- 
tion of the text, "Solution for Pseudosteady States." Since the equations are 
nonlinear, an iterative search procedure employed in the analysis is now 
described. 

For fixed values of the control inputs, angle of attack was increased from 
-10° to 90° in increments of 0.1°. Solutions to equation (11) were obtained 
at each value of a. The real p solutions of equation (11) were used to find 
q, r, and 3 (eqs. (8), (9), and (10)). The resulting values were substituted 
into equation (5) and a residual was computed. If the residual from equa- 
tion (5) changes signs between successive values of a, a common root of equa- 
tions (5) and (11) exists. If the common root was not a root of the determinant 
(eq. (7)), a straight line fit was made between the a values to find an 
approximation to the angle of attack which would yield a zero residual. This 
generally resulted in a residual from equation (5) of less than about 10“^ which 
was considered adequate for determining simultaneous solutions to equations (5) 
and (11). If the residual was greater than 10 ~^, a second straight line fit was 
made between the new value of a and the value associated with the smaller of 
the original two residuals. This procedure was repeated until the residual was 
reduced to less than 10~ 6 . It should be noted that, during the a search, it 
was necessary to track multiple p solution branches to identify which p 
caused the sign change in the residual. 

Having determined the values of p, a, q, r, and 3 which satisfy equa- 
tions (5) and (11), the linearized stability characteristics for perturbations 
about these PSS solutions were calculated. The a search was then continued 
until all PSS solutions were obtained for the specified set of control inputs. 
This search procedure was implemented on both a high-speed computer and a mini- 
computer, and simultaneous solutions of equations (5) and (11) were obtained 
rapidly for any combination of control inputs. 
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TABLE I.- WEIGHT, INERTIA, AND DIMENSIONAL CHARACTERISTICS OF AIRCRAFT 


Weight, W, N 160 968 

Moments of inertia, kg-m 2 : 

I x 35 398 

I y 157 576 

I Z 178 460 

ixz 0 

Wing area, S, m 2 49.2 

Wing span, b,m 11.7 

Mean aerodynamic chord, c,m 4.9 
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TABLE II.- AERODYNAMIC DATA 


[Control derivative units are per degree; all others are per radian] 


a, deg 

C X 

Cx 6 e 

C Z 

Cz «e 

c n(3 

c °p 

Cn r 

Cn <5 a 

Cn 6 r 


Cm 

-10 

-0.0334 

0.00109 

0.721 

-0.00650 

0.1200 

0.0289 

-0.35 

0.00004 

-0.00091 

-2.7 

0.0800 

-5 

-.0469 

.00119 

.360 

-.00650 

.1200 

.0081 

-.35 

.00004 

-.00088 

-2.7 

.0400 

0 

-.0477 

.00102 

-.010 

-.00650 

.1200 

.0015 

-.35 

.00004 

-.00090 

-2.7 

.0100 

5 

-.0320 

.00118 

“.384 

-.00652 

.1200 

.0004 

-.35 

.00004 

-.00081 

-2.7 

-.0200 

10 

-.0173 

.00114 

-.673 

-.00681 

.0970 

.0085 

-.35 

.00007 

-.00068 

-2.7 

-.0380 

15 

-.0177 

.00120 

-.895 

-.00681 

.0745 

.0186 

-.40 

.00015 

-.00064 

-4.0 

-.0590 

20 

-.0442 

.00071 

-1.004 

-.00543 

.0344 

-.0756 

-.60 

.00027 

-.00061 

-5.4 

-.0810 

25 

-.0527 

.00027 

-1.051 

-.00576 

-.1375 

.0079 

-.65 

.00019 

-.00057 

-5.2 

-.1000 

30 

-.0414 

.00009 

-1.195 

-.00514 

-.2177 

.1378 

-.45 

.00019 

-.00044 

-4.5 

-.1330 

35 

-.0430 

-.00058 

-1 .313 

-.00786 

-.2349 

.0506 

-.20 

.00016 

-.00035 

-4.5 

-.1950 

40 

-.0421 

-.00085 

-1.375 

-.00763 

-.1490 

0 

-.10 

.00011 

-.00011 

-5.5 

-.2274 

45 

-.0329 

-.00092 

-1.416 

-.00529 

-.1014 

.0300 

-.05 

.00025 

-.00018 

-6.0 

-.2430 

50 

-.0173 

-.00113 

-1.467 

-.00500 

-.0556 

.0600 

-.04 

.00020 

.00006 

-6.0 

-.2646 

55 

-.0024 

-.00140 

-1.525 

-.00614 

-.0539 

.0596 

-.06 

.00016 

.00001 

-6.0 

-.2871 

60 

.0158 

-.00137 

-1.599 

-.00762 

-.0894 

.0546 

0 

.00017 

.00001 

-5.5 

-.3129 

65 

.0170 

-.00201 

-1.558 

-.00346 

-.1209 

.1151 

.06 

.00007 

.00024 

-3.5 

-.3883 

70 

.0156 

-.00224 

-1.556 

-.00205 

-.1031 

.2147 

.10 

-.00008 

.00010 

-2.0 

-.4304 

75 

.0165 

-.00229 

-1.573 

-.00224 

-.1226 

.5374 

-.30 

.00013 

-.00009 

-2.5 

-.4783 

80 

.0201 

-.00240 

-1.562 

-.00157 

-.2229 

.9108 

-.35 

.00024 

-.00010 

-3.3 

-.5053 

85 

.0254 

-.00249 

-1.557 

-.00133 

-.3747 

.7463 

0 ' 

.00024 

-.00023 

-3.2 

-.5293 

90 | 

.0227 

-.00261 

-1.566 

-.00224 

-.4500 

.6130 

.32 

.00024 

-.00011 

-2.6 

-.5710 



u> 
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TABLE II .- Concluded 


a , deg 

°n>6 e 

( a ) 

% 

c 7 

c 2 

‘r 

C! «, 

u r 

c Yg 

c Y r 

% 


C *6 r 

-10 

- 0.0070 

- 0.0500 

- 0.3010 

- 0.06 

- 0.00070 

0.00018 

- 0.6300 

0.72 

0.0366 

0.00009 

0.00166 

-5 

-.0070 

-.0500 

-.3010 

0 

-.00070 

.00015 

-.6300 

.80 

.0366 

.00009 

.00153 

0 

-.0070 

-.0500 

-.3010 

.03 

-.00070 

.00022 

-.6300 

.82 

.0366 

.00009 

.00159 

5 

-.0070 

-.0917 

-.3014 

.07 

-.00070 

: .00020 

-.6300 

.86 

.0366 

.00009 

.00145 

10 

-.0081 

-.1375 

-.2574 

.15 

-.00064 

.00017 

-.5730 

.84 

.0613 

.00009 

.00116 

15 

-.0088 

-.1260 

-.1022 

.25 

-.00050 

.00011 

i -.5730 

.75 

.1097 

-.00022 

.00124 

20 

-.0072 

-.0570 

-.0991 

.50 

-.00046 

.00006 

-.5160 

.50 

.2894 

-.00032 

.00106 

25 

-.0060 

0 

-.2720 

.62 

-.00027 

.00003 

-.3720 

0 

.3630 

-.00019 

.00132 

30 

-.0055 

.0230 

-.3060 

.48 

-.00020 

.00004 

-.3000 

-.66 

.3158 

-.00015 

.00097 

35 

- 0.0055 

- 0.00670 

0 

-.2526 

.20 

-.00018 

.00002 

-.2900 

-.80 

.1612 

.00018 

.00064 

40 

-.0043 

-.00580 

-.1146 

-.0092 

.07 

.00007 

0 

-.2700 

-.30 

-.5118 

-.00032 

.00049 

45 

-.0043 

-.00485 

-.1146 

-.0380 

.05 

-.00002 

.00007 

-.2600 

-.14 

-.7517 

-.00033 

.00026 

50 

-.0036 

-.00368 

-.1141 

-.0690 

.03 

-.00006 

.00005 

■ -.2500 

-.34 

-.6646 

-.00028 

-.00006 

55 

-.0029 

-.00209 

-.1398 

-.0921 

.01 

-.00001 

0 

-.2400 

-.38 

-.6069 

0 

-.00020 

60 

0 

-.00125 

-.1501 

-.0965 

0 

0 

.00003 

-.2200 

-.38 

-.5267 

-.00014 

.00001 

65 

0 

-.00311 

-.1455 

-.0805 

0 

.00003 

.00002 

.0499 

-.38 

-.7643 

-.00022 

.00016 

70 

0 

-.00329 

-.1553 

-.1381 

.01 

.00007 

-.00001 

.1720 

-.27 

- 1.0630 

-.00003 

.00019 

75 

0 

-.00408 

-.1719 

-.1991 

.02 

.00004 

.00002 

.1720 

-.60 

-.9094 

-.00010 

.00013 

80 

0 

-.00335 

-.1845 

-.1815 

.02 

-.00007 

0 

-.0880 

-.90 

-.2661 

-.00016 

.00020 

85 

0 

-.00349 

-.1868 

-.2040 

.02 

-.00008 

.00009 

-.2740 

-.90 

-.6982 

-.00013 

.00014 

90 

0 

-.00350 

-.1900 

-.1990 

0 

-.00005 

.00004 

-.3500 

- 1 .00 

- 1.1560 

-.00010 

.00020 


a For a = 35 ° to 90 °, the left column of values is for 6 e ^ 0 f and the right column for 6 e < 0 








TABLE IV.- STABILITY CHARACTERISTICS OF 


ONE-G TRIM SOLUTIONS 


Curve 

6 a , deg 

Characteristic roots 

IP 

0 

-0.69 

-0.31 ± 1 . 50i 

-0.24 ± 2 . 1 i 


-15 

-.62 

-.36 ± 1 . 40i 

-.25 ± 2. 7i 


-30 

-.81 

-.25 ± 1 . 70i 

-.11 ± 3. Oi 

2N 

0 

2.00 

-1 .38 ± 0 . 56i 

-.55 ± 2. 9i 


-15 

1.60 

-.68 ± 1 .60i 

-.89 ± 3 . Oi 


-30 

.53 

-.57 ± 1 . 60i 

-.41 ± 3 . 2i 

2P 

0 

2.00 

-1.38 ± 0 . 56i 

-.55 ± 2. 9i 


-15 

2.20 

-1.71 ± 0.48i 

-.19 ± 3.6i 


-25 

2.30 

-1.78 ± 1 . 30i 

-.13 ± 3 . 7i 

3P 

0 

-2.20 

1.58 ± 2. 50i 

-1.00 ± 3.4i 


-15 

-2.30 

1 .42 ± 2. 20i 

-.88 ± 3 . 2i 


-25 

-2.70 

1.33 ± 1 .70i 

-.69 ± 3 . Oi 

3N 

0 

-2.20 

1 .58 ± 2 . 50i 

-1.00 ± 3.4i 


-15 

-2.30 

1.71 ± 2 . 70i 

-1.10 ± 3 . 6i 


-30 

-2.40 

1.82 ± 2. 80i 

-1.10 ± 3. 7i 

4N 

0 

.40 

a_ . 91 ± 3 . OOi 

.25 ± 3 . 6i 


-15 

1.10 

-1.30 ± 2 . 90i 

.41 ± 3. 5i 


-30 

1.40 

-1.50 ± 2. 90i 

.42 ± 3 . 6i 

4P 

0 

.40 

-.91 ± 0 . 30i 

.25 ± 3. 6i 


-15 

.30 

-.54 ± 3 . 20i 

-.10 ± 4. 1 i 


-30 

1 .00 

-.72 ± 3 . 90i 

-.17 ± 4. 5i 

5N 

0 

-1.90 

-.33 ± 5 . 70i 

.90 ± 3 . 5i 


-15 

-1.80 

-.32 ± 5 . 70i 

.83 ± 3. 4i 


-30 

-1 .60 

-.32 ± 5 . 70i 

.74 ± 3.3i 

5P 

0 

-1.90 

-.33 ± 5 . 70i 

.90 ± 3. 5i 


-15 

-2.10 

-.33 ± 5 . 70i 

.97 ± 3. 5i 


-30 

-2.20 

-.33 ± 5 . 70i 

1.04 ± 3 . 6i 

6N 

0 

2.80 

-.41 ± 7 . 70i 

1.43 ± 3. 1 i 


-15 

2.70 

-.40 ± 7 . 40i 

1.43 ± 3.1 i 


-30 

2.70 

-.39 ± 7 . 30i 

1.44 ± 3. 1 i 

6P 

0 

2.80 

-.41 ± 7 . 70i 

1.43 ± 3. 1 i 


-15 

2.80 

-.41 ± 8. OOi 

1.43 ± 3.1 i 


-30 

2.80 

-.42 ± 8 . 20i 

1.42 ± 3.1 i 


TABLE V.- STABILITY CHARACTERISTICS OF BASIC AND ASSOCIATED SOLUTIONS 


A6 e , deg 

<5 a r de 9 

Characteristic roots 

8 

0 

-0.92 

-0.130 

± 1 . 70i 

-0.36 ± 2.11 


-4 

-.92 

-.190 

± 1 . 40i 

-.29 ± 2 . 4i 


-8 

-.94 

-.210 

± l.OOi 

-.27 ± 2 . 7i 


0 

0.63 

-1.100 

± 0 . 42i 

-0.10 ± 3. 9i 


-4 

.57 

-1.000 

± . 38i 

-.18 ± 3 . 9i 


-8 

.46 

-.960 

± .311 

-.20 ± 3 . 7i 


0 

-0.86 

-0.200 

± 1 .601 

-0.43 ± 3 . 9i 


-15 

-.97 

-.190 

± 1 .30i 

-.36 ± 4 . 1 i 


-30 

-.84 

-.310 

± 1.201 

-.30 ± 4.11 

6 

0 

-0.90 

-0.150 

± 1 . 70i 

-0.36 ± 2.11 


-5 

-.90 

-.220 

± 1 .401 

-.28 ± 2.51 


-10 

-.90 

-.240 

± .801 

-.24 ± 3.01 


-5 

0.54 

-0.21 

-1.80 

-2.40 ± 3. 9i 


-10 

.20 

-.65 

-.97 

-2.10 ± 3 . 8i 


-5 

-1 .40 

-0.100 

+ 0 . 48i 

o 

+i 

r- 

CM 

o 


-15 

-.84 

-.340 

± . 801 

-.30 ± 4.11 


-30 

-.82 

-.310 

± . 97i 

-.28 ± 4 . 8i 

4 

0 

-0.85 

-0.170 

± 1 . 80i 

-0.36 ± 2.11 


-15 

-.65 

-.370 

± .601 

-.27 ± 3. 4i 


-25 

-.73 

-.320 

± .341 

-.27 ± 3. 8i 


-27 

-.85 

-.270 

± .021 

-.24 ± 4 . 3i 


-30 

-.87 

-.270 

± .651 

-.24 ± 4 . 7i 

2 

0 

-0.78 

-0.200 

± 1.901 

-0.36 ± 2.11 


-15 

-.66 

-.370 

± 1 .101 

-.30 ± 3.21 


-30 

-. 95 

-.300 

± 1.101 

-.18 ± 3.31 

0 

0 

-0.69 

-0.310 

± 1 .50i 

-0.24 ± 2.11 


-5 

-.63 

-.340 

± 1 .301 

-.23 ± 2 . 3i 


-15 

-.62 

-.360 

± 1.401 

-.25 ± 2 . 7i 


-30 

-.81 

-.250 

± 1 . 70i 

-.11 ± 3.01 

-4 

0 

-0.25 

-0.280 

± 2 . 20i 

-0.18 ± 2 . 5i 


-15 

-.45 

-.280 

± 2.101 

-.14 ± 2.81 


-30 

-.75 

-.200 

± 1.801 

-.10 ± 3 . 2i 

-8 

0 

-0.35 

-0.170 

± 2 . OOi 

-0.36 ± 2.41 


-15 

-.37 

-.320 

± 1.701 

-.30 ± 2.81 


-30 

-.42 

-.230 

± 1.301 

-.28 ± 3 . 5i 

-12 

0 

-1.70 

0.330 

± 0.521 

-0.31 ± 2.31 


-1 

-1 .80 

.360 

± . 41 i 

-.32 ± 2 . 3i 


0 

-2.50 

1.6 

-0.22 

-0.43 ± 2.61 


-1 

-2.30 

1.4 

-.26 

-.35 ± 2 . 4i 



w 


TABLE VI.- COMPARISON OF PSS SOLUTIONS WITH EQUILIBRIUM SPIN CONDITIONS (A6 e = 6 a = 6 r = 0°) 


Method 

a, 

deg 

P/ 

deg/sec 

<3/ 

deg/sec 

r, 

deg/sec 

6, 

deg 

V, 

m/sec 

Characteristic roots 

Reference 1 1 

35.7 

■EH 



-0.90 

187 

-2.00 

-0.10 

-0.05 ± 1 . 6i 

1.40 ± 2. Oi 

-0.80 ± 2. 5i 

PSS 

37.5 


mEm 


-1.20 

266 

-2.20 



1.60 + 2.5i 

-1.00 ± 3.4i 

PSS (D = W) 

37.6 

ID 

ifsfD 


-.86 

181 

-1.50 



1.20 ± 1 . 7i 

-.80 + 2. 1 i 

PSS (D = w. 

37.4 

m 



-1 .60 

181 

-1.70 



1.30 ± 1.7i 

-.80 ± 2.1 i 

Y = -90°) 












Reference 11 

60.4 

49.3 

3.40 

86.5 

0.50 

142 

-0.13 

0.04 

-0.07 ± 1 . 7i 

0.20 ± 2. 6i 

-0.40 ± 2 . 4i 

PSS 

49.3 

100.2 

1.30 

116.3 

.40 

266 

.40 



.30 ± 3. 6i 

-.90 ± 3 . Oi 

PSS (D = W) 

46.6 

62.1 

-.30 

65.4 

0 

162 

.39 



.30 ± 2. Oi 

-.70 ± 1.6i 

PSS (D = W, 

52.2 

56.7 

1.10 

72.8 

-1.00 

153 

.18 



.04 ± 2. li 

-.40 ± 1 .5i 

Y = -90°) 












Reference 1 1 

73.1 

40.9 

1.60* 

133.5 

0.10 

136 

-0.42 

0.15 

-0.07 ± 2.4i 

0.30 ± 2. 9i 

-0.20 ± 2.8i 

PSS 

73.3 

79.2 

-3.70 

263.4 

-.90 

266 

-1.90 



.90 ± 3. 5i 

-.30 ± 5 . 7i 

PSS (D = W) 

73.6 

40.2 

-3.80 

135.5 

-.40 

136 

-1.10 



.50 ± 1.6i 

-.20 ± 2 . 9i 

PSS (D = W, 

73.3 

40.6 

-1.20 

134.4 

-1.00 

136 

-1.10 



.50 ± 1 . 6i 

-.20 ± 2. 9i 

Y = -90°) 












Reference 11 

82.1 

28.9 

1.50 

206.0 

0.23 

134 

0.27 

-0.15 

-0.06 ± 3.6i 

-0.40 ± 3.4i 

-0.01 ± 3. 9i 

PSS 

83.7 

51.6 

-.80 

460.7 

.03 

266 

2.80 



-1.40 ± 3. 1 i 

-.40 ± 7.7i 

PSS (D = W) 

83.9 

25.5 

-1.60 

235.3 

.50 

133 

' 1.40 



-.70 ± 1 . 5i 

-.20 ± 3.9i 

PSS (D = W, 

83.7 

26.1 

.01 

230.4 

-.13 

134 

1.40 



-.70 ± 1 . 5i 

-.20 ± 3 . 9i 

Y = -90°) 


















P Positive roll-rate solutions 

N Negative roll-rate solutions Stable 


Curves numbered 1 to 6 in order Unstable oscillation 

of increasing angle of attack Divergent 




(a) Angle of attack. 

Figure 2.- PSS solution for one-g trim condition A6 e = o°. 





30 20 10 0 -10 -20 


6 a , deg 

(c) Yaw rate. 
Figure 2.- Continued. 
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p, deg/sec 



Figure 3.- Basic and associated PSS solutions with various elevator 
and aileron inputs for ot < 25°. 







q, deg/sec 



Stable 

Unstable .oscil lation 

Divergent 


-401 1 1 L - • — J - 1 -- 
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6 a > deg 

(d) Pitch rate. 


Figure 3.- Continued. 





200 


(D to (6) indicate PSS predictions where numbers 
refer to solution types shown on figure 2 



a , deg 


(a) Roll rate; 


6 = ± 30 ° 

a 

Q 

© 

30 ° 


70 


80 


Figure 4.- Roll-rate and angular-velocity PSS solutions as a function of angle of attack, 
^a.min t ' ie ^a at wtl ich solutions 2 and 3 join (see fig. 2). 
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PSS prediction 
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Figure 8.- Calculated responses initiated from pitch-down condition 
for S-, = -15°. A 6 = = 8 °: a = -3.8°. 
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Time, sec 

(b) Dominant components of r. 



1 2 
Time, sec 

(c) Dominant components of p 


Figure 11.- Concluded. 



H] to GO indi cate PSS predictions where 
numbers refer to solution types on figure 2 



Figure 12.- Phase plots for pitch-up maneuver. AS 


= - 12 °; 




500 


\z\ to [6] indicate PSS predictions where numbers 
refer to solution types shown on figure 2 


O Unstable oscillation 
□ Divergent 


© 


400 



(c) w and a phase plot. 
Figure 12.- Concluded. 
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